MATH 18.152 COURSE NOTES - CLASS MEETING # 13 



18.152 Introduction to PDEs, Fall 2011 Professor: Jared Speck 

Class Meeting jj^ 13: Geometric Energy Estimates 



1. Dm, THE ENERGY-MOMENTUM TENSOR, AND COMPATIBLE CURRENTS 

The following shorthand notation is often used for the "linear wave operator associated to m:" 

(1.0.1) u^^i'im-'r^d^dp. 

Using this notation, the wave equation -d"^ + A0 = can be expressed as 

(1.0.2) □,n0 = O. 

We now introduce a very important object called the energy-momentum tensor. As we will see. 



it encodes some very important conservation laws associated to solutions of (1.0.2) 



Definition 1.0.1. The energy- momentum tensor associated to equation (1.0.2) is 



(1.0.3) T^u^= d^<pd,<P-]^m^,{m-^T^d^<pdp4). 

Later in the course, we will hopefully have time to motivate its derivation in the larger context 
of variational methods. For now, we will simply study/use its useful properties. 
Note that T^y is symmetric: 

(1.0.4) T^u = n^. 

In your homework, you will prove the following very important positivity property of T, which is 
called the dominant energy condition. 

Lemma 1.0.1 (Dominant Energy Condition for T^^). 
(1.0.5) 

T{x,Y) =-^T„^x"r^ > are both timelike and future-directed or timelike and past- directed. 



Since causal vectors are the limit of timelike vectors, we have the following consequence of (1.0.5): 
(1.0.6) 

T(X, Y) = T^pX'^Y^ > if X, F are future-directed and causal or past-directed and causal. 
As before, we can raise the indices of T : 

(1.0.7) T'^^ = {m-^y{m-^y^To,p. 
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A very special case of Lemma 1.0.1 is the following, which corresponds to X^^ - - 6q 
(1, 0, 0, 0) in the lemma: 

(1-0.8) Too = T'' = lt(d,<pr = l\Vt,.<P\'. 

^ /t=0 ^ 



The derivation of (1.0.8) is a simple computation that you should do for yourself. Note that Too is 



positive definite in all of the derivatives of 0. This fact will play an important role in Theorem |2.1 
below. 

The next lemma shows that T^" is divergence-free whenever verifies the wave equation. This fact 
is intimately connected to the derivation of conservation laws, which are fundamental ingredients 
in the study of hyperbolic PDEs. 

Lemma 1.0.2 (The divergence of T^^^). Let T^^, he the energy-momentum tensor defined in 



(1.0.3). Then 



1.0.9) d;r^- = {u^4>)ijn-^Td^4>. 



In particular, if (p is a solution to (1.0.2), then 



(1.0.10) df,T^'' = 0. 

Proof. The proof is a computation that uses the symmetry property {m'^)^^'^ = {^n^^Y^^ and the 
fact that we are allowed to interchange the order of partial derivatives (if is sufficiently smooth) : 

(1.0.11) d^T^''' = d^i^{rn-^Y''{rn-^y^do,(pdp(p - ]^{rn-^Y'' {rn^Y^dc 

= {n^<f)){m''Y^dp<j)+{m-'Y''{m-'Y''ida<P)df.dp<l> 

= in^<j))im-'Y^d^<P, 
where the last three terms have canceled each other. 

□ 

As we will soon see, the energy-momentum tensor provides an amazingly convenient way of 
bookkeeping in the divergence theorem. However, in order to apply the divergence theorem, we 
need to find a useful vectorfield to take the divergence of. By useful, we mean a vectorfield that can 
be used to control a solution to the wave equation. One way of constructing a useful vectorfield is 
to start with an auxiliary vectorfield X and then to contract it with the energy momentum tensor 
to form a new vectorfield J. The next definition shows how to do this. 

Definition 1.0.2. Given any vectorfield X, we associate to it the following compatible current, 
which is itself a vectorfield: 

(1.0.12) Wj^=^T^"X„. 
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So which vectors X are the useful ones? It turns out that the answer is causal vectors. This fact is 



closely connected to the dominant energy condition (|1.0.5|). This will become more clear in our proof 
of theorem 



2.1 



below; Note that by Lemma 



1.0.1 



To.pX'^YP =^ T(X, F) > if 



X, y are both timelike and future-directed (i.e., X°,y° > 0) or past-directed (i.e., X'',y° < 0). 

In order to apply the divergence theorem to (^)J^, we of course need to know its divergence. We 
carry out this computation in the next corollary. 



Corollary 1.0.3. Using (1.0.4) and (1.0.10), we have that 



(1.0.13) 

where 



(1.0.14) 

is called the deformation tensor of X. 



defl. 



"r^t^u = -(d^,X^ + dyX^ 



We now state a version of the divergence theorem that is tailored to our study of the linear wave 
equation. 

Theorem 1.1 (Divergence Theorem). Let cj) he a solution to the linear wave equati on Om4 ' = 0- 
Let X be any vectorfield, and let ^^^J be the compatible current define in Definition 1.0.2. Let 
Vt c ]Ri+" be a domain with boundary dVL. Then the following integral identity holds: 



1.0.15) 



JdQ. Jn ^ ' 

2. Energy Estimates and Uniqueness 



x. 



We will now use the results of the previous section to derive some extremely important energy 
estimates for solutions to = 0. The results we derive are a geometry version of integration by 
parts + the divergence theorem. They could alternatively be derived by multiplying both sides of 
the wave equation by a suitable quantity and then integrating by parts over a suitable hypersurfaces, 
but there is a substantial gain in geometric insight that accompanies our use of compatible currents. 

Theorem 2.1 (Energy estimates in a cone). Let (f)(t,x) be a solution to the 1+n dimensional 
global Cauchy problem for the linear wave equation 



(2.0.16) 
(2.0.17) 
(2.0.18) 



□r 



0, 



^iO,x) = f(x), 
dt(j){0,x) =g{x), 



X € 



X € 



Let Re [0, cx)], let X be the past- directed timelike vector defined by X>^ = -6q, and let ^^^J^^[(f){t,y)] 
be the compatible current (1.0.12) associated to X. Note that by (1.0.8), (^^J^[0(t, y)] = |Vt,j,0(t, y)p - 



E)!=o(^M0)^ = (dt4>y + TJi=i{di(j)Y . Define the square of the energy by 
(2.0.19) E\m)'= [ N,^''^nct>{t,y)]ry = \ f \Vt,yHt,y)\' d^V, 
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where = 6^ (and therefore iV^ = -6q) is the past-pointing unit normal covector to {t} x Bn^tip) 
M^, and Br{p) c denotes the solid Euclidean hall of radius R centered at p. Then 



(2.0.20) 



i5;[0](t)<£;[0](o). 



Proof. The goal is to apply Theorem 1.1 to the solid truncated backwards light cone Ct^p-R {(r, y) e 
[0, oo) X M" I \y - p\ < R - t} and to make use of the dominant energy condition. It is easy to see 

that dCt^p-R = B u M.t,p-R u T, where B {0} x Br{p) is the fiat base of the truncated cone, T 

{t] X BR^t{p) is the fiat top of the truncated cone, and M.t,p-R {(t, y) e [0, oo) xM" | \y-p\ = R-t} 
is the mantle of the truncated cone. 



By Theorem 1.1, we have that 



(2.0.21) 
where 



E[<pm-E[<pm + F[<p]= / dj J^[<PiT,y)])dTd-y 

■JCt.v 



(2.0.22) 



is the "flux" associated to M.t p-R- Since 



'Mt,p-R 

solves the wave equation (2.0.16), and since ("^^vr 



the identity (1.0.13) implies that the right-hand side of (2.0.21) is 0. Therefore, 



(2.0.23) E[(j)](t) - E[(j)]{0) + F[(l)] = 0. 

We claim that F[(j)] > 0. The energy inequality (2.0.20) would then follow from (2.0.23). They key 
observation for showing that F[(j)] > is the following. Along the mantle Mt,p-R, it is easy to see 
(draw the picture!) that = Llu-, where L is a past-directed null vector. Therefore, the integrand 



in (2.0.22) is equal to TafsX""!/, and since X is a past-directed timelike vector, the dominant energy 



condition (1.0.6) implies that T^isX^'lJ > 0. Therefore, F[(j)] > as desired. 



□ 



1 and 02 to (2.0.16) have the 
Then the two solutions agree on the 



Theorem |2.1| can easily be used to prove the following local uniqueness result for solutions to the 
linear wave equation. 

Corollary 2.0.4 (Uniqueness). Suppose that two solutions 
same initial data on Br(p) c{(r, y) | r = 0} 

light cone" Cp.R=^ {{T,y) \ < r < R,0 <\y - p\ < R- t}. 
Proof. Define 'ijj'^=^ (pi 



"solid backwards 



-02- Then ip verifies ( |2.0.16[ ) and furthermore, i?['0](O) = 0. Thus, by Theorem 
2.1 , E[ilj]{t) = for <t < R. Therefore, from the definition of E[il)]{t), it follows that Vr,y'0'(T, y) = 



for (r, y) e Cp.R. Thus, by elementary analysis, ip is constant in Cp-^R. But -0(0, x) = for (0, x) e Cp-R. 
Thus, ip(T,y) = for all points (t, ?/) € Cp-R. □ 



Corollary 2.0.4 is one illustration of the finite speed of propagation property associated to the 
linear wave equation. Another way to think about it is the following. Suppose you alter the initial 
conditions outside of Br{p), but not on Br(p) itself. Then this alteration has no effect whatsoever 
on the behavior of the solution in the spacetime region Cp-R Think about this claim yourself; it 
follows easily from the Corollary! 
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3. Developments, Domain of Dependence, and Range of Influence 



We will now develop a language for discussing the finite speed of propagation properties of the 
linear wave equation in more detail. If we had more time in this course, we could adopt a more 
geometric point of view that would apply to many other hyperbolic PDEs. This would involve 
fleshing out our discussion of Lorentzian geometry, and also developing a generalized version of 
geometry that applies to a large class of PDEs. 

Warning 3.0.1. Some people permute or even severely alter the following definitions, which can 
be very confusing. The definitions below therefore indicate some of my biases. 

Definition 3.0.3 (Development). Let S c {(t,x) | t = 0} be a set. Assume that that we know the 
initial data 0(0, - f{x), dt(f){0,x) = g{x) for the wave equation (1.0.2), but only for x e S. Then a 
future development Q of S is defined to be a "future" region of spacetime Q c ]Ri+" n {{t,x) | t > 0} 



on which the solution (f){t,x) to (1.0.2) is uniquely determined by the initial data on S. A past 
development 'D'(S) can be analogously defined (replace t>0 with t < in the previous definition). 



Example 3.0.1. If Bji(p) and Cp-R are as in Corollary 2.0.4, then Cp.ji is a development of B^^p). 
You can imagine that the solution knows how to "develop" in Cp-^R from the initial conditions on its 
subset Br{p). 

Definition 3.0.4 (Maximal development). The maximal future development of S, which we 
denote by 'D^{S), is defined to be the union of all future developments of S. The maximal past 
development 'D'{S) can be analogously defined. The maximal development of 5* is defined to be 
V+(S)uV-(S). 



Example 3.0.2. Consider the plane P {{t,x^ 







x^ 



0}. Then using techniques from 



a more advanced course, one could show that 'D(P) = P for the wave equation (1.0.2). That is, 
knowing the conditions of a solution cj) along P is not enough information to determine the solution 
anywhere else. This is closely connected to the fact that all smooth curves in P have tangent vectors 
that are timelike relative to the Minkowski metric. 

Definition 3.0.5 (Domain of dependence). Let Vt c ]Ri+". Assume that is a solution to the 
wave equation (1.0.2) in f2. A domain of dependence for is a set S such that is completely 
determined on Vt from only the data 01^ and Vt,x'P\s- 

Remark 3.0.1. For general nonlinear hyperbolic PDEs, domains of dependence depend both on 
Vt and the solution itself. However, for the linear wave equation, domains of dependence do 
not depend on the solution. Roughly speaking, this is because the "geometry of the solution" is 
predetermined by the Minkowski metric m. 

Example 3.0.3. In 1 + 1 dimensions, a domain of dependence for the spacetime point (t, x) (for the 



wave equation (1.0.2)) is the "initial data" interval {0} x \_x-t,x + i\. Another domain of dependence 
for this point is the interval {t/2} x [x - t/2, x + t/2]. A trivial example is that (t, x) is a domain of 
dependence for itself. 

Example 3.0.4. In 1+3 dimensions, a domain of dependence for the positive t axis {(t, x^, x^, x^) | x^ : 
2 - _ g^^ > g| ^fQj^ wave equation (1.0.2)) is all of "space:" {(t, x^, x^, x^) | t = 0}. Any subset 



X 



X"" 



of space is not a domain of dependence for the positive t axis. 



The next definition is complementary to the notion of domain of dependence. 
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Definition 3.0.6 (Range of infiuence). Assume that is a solution to the wave equation (1.0.2) 
in The range of influence TZ for a set S c ]Ri+" is the set of all points € ]R^+" such that 

(p{t,x) is affected by the initial data 0|s and Vt,a;0|5- 

Example 3.0.5. In 1 + 1 dimensions, the (future) range of infiuence (for t > 0) of the interval 
5 = {0}x[-l,l] isn = {{t,x) I -t-l<x<t+l}. 

Example 3.0.6. In 1 + 1 dimensions, the (future) range of influence (for t > 0) of the t axis 
S = {(t,0) \t>0} is7^ = {(t,x) I t>0}. 

Example 3.0.7. In 1 + 3 dimensions, the (future) range of influence (for t > 0) of S" = {0} x dBi(0) 

isn = {(t,x) I <t < = 1 -t} u{(t,x) \0<t<t,\x\ = l+t} where |x| ^= y/ (x^y + (x^y + (x^y . 
This is a consequence of the Sharp Huygens' Principle. 
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